arXiv:1501.03652vl [cond-mat.mes-hall] 15 Jan 2015 


Aharonov Bohm effect in 2D topological insulator. 
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We present magnetotransport measurements in HgTe quantum well with inverted band structure, 
which expected to be a two-dimensional topological insulator having the bulk gap with helical gapless 
states at the edge. The negative magnetoresistance is observed in the local and nonlocal resistance 
configuration followed by the periodic oscillations damping with magnetic field. We attribute such 
behaviour to Aharonov-Bohm effect due to magnetic flux through the charge carrier puddles coupled 
to the helical edge states. The characteristic size of these puddles is about 100 nm. 


I. INTRODUCTION 

The investigation of the quantum interference phenom¬ 
ena, such as Aharonov-Bohm (AB) oscillations and weak 
localizations, provides an important information about 
fundamental properties of various electronic systemai. 
Recently new class of materials with interesting prop¬ 
erties have emerged, called topological insulators (TI), 
which are insulating in the bulk and characterized by the 
existence of robust gapless excitations at their surface^^— . 
Manifestation of AB oscillations in topological insulators, 
can be utilized to probe the surface nature in TI and its 
properties. The time-reversal-symmetric 2D topological 
insulator is induced by a strong spin-orbit interaction^AS 
and characterized by edge modes with opposite spins 
propagating in opposite directions. The 2D TI have 
been realized in HgTe quantum wells with inverted band 
structureiid^. Aharonov Bohm oscillations is intimately 
related to the weak localization (WL) corrections. For 
example, WL phenomena is discussed in terms of the col¬ 
lective action of the magnetic flux through the random 
loops of the pair of time reversed trajectories^. The con¬ 
structive interference of these path in return point is sup¬ 
pressed by AB flux, which lead to the resistance increase 
and negative magnetoresistancei^. The spin orbit cou¬ 
pling strongly modifies the quantum corrections the con¬ 
ductivity, it leads to a destructive interference between 
clockwise and counter-clockwise trajectories and changes 
the WL to weak anti-localization behavior (WAL). 

In 2D topological system WL is quite different from 
conventional 2D metals and strongly affected by the 
Dirac spectrum of massive fermion a^^d^ with the mass 
proportional to the band gap. The electronic spectrum 
is shown in Figure 1. The band structure has a sim¬ 
ple parabolic form near the bottom of the band (marked 
by Ep in Figure 1 and linear Dirac-like dispersion re¬ 
lation at high Fermi energy Ep. When the Fermi en¬ 
ergy lies in the bulk gap near the charge neutrality point 
(CNF), spectrum described by pair of the helical edge 
states (Figure 1). More over, the edge state must have 
linear Dirac like dispersion. The localization corrections 
are strongly governed by such remarkable property of the 
spectrum. When the Fermi energy is small, the effec¬ 


tive spin-orbit coupling is weak, and one can expect the 
conventional WL behaviour—. When the Fermi energy 
becomes larger than the gap width, the energy disper¬ 
sion is linear, and theory describes crossover from WL to 
WAL behaviour— A2. The magnetoresistance of 2D TI 
with a dominant edge state contribution is described by 
two mechanisms: first scenario relies to the frequent de¬ 
viations of the edge electrons into the disordered AB flux 
threaded 2D bulbil; and second scenario to the localiza¬ 
tion of the helical edge sates due to the collective action 
of the random magnetic flux through the loops, naturally 
formed by rough edges of the sampleiS. Both models 
predict quasi-linear positive magnetoresistance MR. Ac¬ 
cording to this theoretical predictions positive MR was 
observed in 8 nm HgTe quantum wells^Si. 
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FIG. 1: (Color online)Schematic of the energy spectrum of 
2D TI near k=0. Solid lines -spectrum for the bulk electrons, 
dashes- helical edge state spectrum. Positions of the Fermi 
energy in the bulk gap and near the bottom of the conductive 
bands are indicated. 

Recently interaction of the helical states with multi¬ 
ple puddles of charge carriers formed by fluctuations in 
the donor density has been considered in 2D topological 
insulators2i. This leads to significant inelastic backscat- 
tering due to the dwelling of carriers in the puddles and 
to a large resistance, which depends weakly on temper- 
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ature. Moreover, carriers in the puddles describes loops 
trajectories, which are sensitive to the magnetic field flux, 
therefore, one could expect that the magnetoconductiv¬ 
ity is mainly contributed by the quantum interference of 
such electrons. 

In the present paper we investigate transport prop¬ 
erties of the HgTe quantum wells with the width d of 
8-8.3 nm. Slightly above the the charge neutrality point, 
when the system is expected to be two-dimensional topo¬ 
logical insulator, we observe the negative magnetoresis¬ 
tance measured in the local and nonlocal configurations, 
followed by periodic oscillations damping with magnetic 
field. When the applied field is larger than 1 T, the MR 
becomes positive. From the sign of the magnetoresistance 
we specified two contributions: electrons from edge states 
and bulk electrons, which are localized in the metallic 
puddles formed by potential fluctuations. 


II. EXPERIMENT 


The Cdo.e 5 Hgo. 35 Te/HgTe/CdQ,e 5 Hgo. 35 Te quantum 
wells with (013) surface orientations and a width d of 8- 
8.3 nm were prepared by molecular beam epitaxy. A de¬ 
tailed description of the sample structure has been given 
i n2^i23 . Device A is six-probe Hall bar, while device B is 
designed for multiterminal measurements. The device A 
was fabricated with a lithographic length 6/rm and width 
5^m (Figure 2, top panel). The device B consists of 
three Ag,m wide consecutive segments of different length 
(2,8,32^m), and 7 voltage probes. Device C (figure 5) 
is a structure with large gate area for identifying nonlo¬ 
cal transport over macroscopic distances^^. The lengths 
of the edge states are determined by the perimeter of 
the sample part covered by metallic gate (mostly side 
branches) rather than by the length of the bar itself. The 
ohmic contacts to the two-dimensional gas were formed 
by in-burning of indium. To prepare the gate, a dielectric 
layer containing 100 nm Si 02 and 200 nm Si^Ni^ was 
first grown on the structure using the plasmochemical 
method. Then, the TiAu gate with sizes of 18 x lO^m^ 
was deposited. The ungated HgTe well was initially n- 
doped with density Ug = 1.8 x Several devices 

with the same configuration have been studied. The den¬ 
sity variation with gate voltage was 1.09 x . 

The magnetotransport measurements in the structures 
described were performed in the temperature range 1.4- 
25 K and in magnetic fields up to 12 T using a stan¬ 
dard four point circuit with a 3-13 Hz ac current of 0.1- 
10 nA through the sample, which is sufficiently low to 
avoid overheating effects. The density of the carriers in 
the HgTe quantum wells can be electrically manipulated 
with local gate voltage Vg. The typical dependence of 
the four-terminal resistance of two of the representative 
samples A and B as a function of Vg is shown in Figure 
2. The resistance i?i 4,23 = 7?/=i,4;iz=2,3 of the sample A 
and resistances for sample B, measured by various volt¬ 
age probes in a zero magnetic field reveal a sharp peaks, 
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FIG. 2: (Color online) Color online) (a) a) The local and 
nonlocal resistances A as a function of the gate voltage at zero 
magnetic field measured by various voltage probes for samples 
A (a) and B (b), T=4.2 K. Top panel-shows schematic view 
of the samples. The perimeter of the gate is shown by blue 
rectangle. 


shown in figure 2 a and b, when the gate voltage induces 
an additional charge density, altering the quantum wells 
from an n-type conductor to a p-type conductor via a 2D 
TI state. It has been show n^^d^ that the 4-probe resis¬ 
tance in an HgTe/CdTe micrometer-sized ballistic Hall 
bar demonstrated a quantized plateaux Ru ^23 — hl2e^. 
It is expected that the scattering between the helical edge 
states in the topological insulator is unaffected by the 
presence of a weak disorderi^i^ii. Note, however, that the 
resistance of samples longer than 1/im might be much 
higher than due to the presence of the electron spin 

flip backscattering on each boundary. Mechanism of the 
back scattering is not clear and appealing task for theo¬ 
reticians and a matter of ongoing debat o^^d^d^ . The Hall 
effect reverses its sign and Rxy ~ 0 (not sown) when lon¬ 
gitudinal resistance Rxx approaches its maximum value, 
which can be identified as the charge neutrality point 
(CNP). These behaviour is similar to those described in 
graphene^l. An unambiguous way to prove the presence 
of edge state transport mechanism in 2D TI with strong 
backscattering on the boundary are the nonlocal electri¬ 
cal measurements. The nonlocal response always exists 
because of the presence of the two counter-propagating 
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edge states, which flows sideways and may reach any con¬ 
tacts in the device^ii^^. Figures 2 a and b show the 
nonlocal resistances corresponding to the different con¬ 
figurations. For example, nonlocal resistance in figure 2a 
corresponds to the contact configuration, when the cur¬ 
rent flows between contacts 2 and 6 and the voltage is 
measured between contacts 3 and 5. One can see that 
the nonlocal resistance near CNF has a peak of a com¬ 
parable amplitude, though less wide, and approximately 
in the same position as the local resistance. Outside of 
the peak the nonlocal resistance is negligibly small. The 
apparent residual nonlocal resistance in figure 2b above 
Vg — VcNP > and below Vg — Vcnp < —11^ is is 
related to the noise in the full linear scale measurements. 
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FIG. 3: (Color online) Color online) The local resistance R 
for two samples with multiterminal configuration (devices B) 
as a function of the magnetic field for the different values 
of the gate voltages, when the Fermi energy passes through 
CNP from the hole to electronic side of the resistance peak, 
(a) Vg — Vcnp (V): -0.15 (cyan), -0.02 (blue), 0.1 (green), 
0.35 (red), 0.47 (black). (h)Vg — Vcnp (V): -0.18 (cyan), 0.05 
(green), 0.25 ( blue), 0.37 (red), 0.45 (black). 

It is expected that a magnetic field perpendicular to 
the quantum well breaks time reversal symmetry (TRS) 
and thereby enables elastic scattering between counter- 
propagating chiral edge modes. However, a number 
of the different theoretical models has previously been 
propose d^^d^i^^i^^i^o with substantially different physical 
scenarios. In our previous study^S we have observed a lin¬ 


ear negative magnetoconductance in HgTe-based quan¬ 
tum wells in the 2D TI regime near CNP, when the edge 
state transport prevails. Our observation agrees with the 
model^ which describes the effects of WAL (see discus¬ 
sion section below for more details). The model predicts 
almost linear positive magnetoresistance = A|i?|, 
where parameter A strongly depends on the disorder 
strength W in comparison with the energy gap Eg : the 
fluctuations W > Eg result in a large B-slope corre¬ 
sponding to a strong disorder regime and fluctuations 
W < Eg lead a small B-slope. Figure 3 shows the evolu¬ 
tion of the resistance Rxx with magnetic field and den¬ 
sity, when the chemical potential crosses the bulk gap. 
The magnetoresistance (MR) demonstrates a striking V- 
shape dependence in magnetic fields below IT near CNP, 
which confirms our previous observations. Note however, 
that the V-shaped magnetoresistance is strongly trans¬ 
formed, when Fermi energy moves away from CNP to 
the electronic side of the resistance peak. In this region 
MR inverses the sign near zero magnetic field and shows 
triangular-shaped peak accompanied by two satellite fea¬ 
tures or damped oscillations. The width of the negative 
magnetoresistance spike and the period of the oscillations 
are slightly varied from sample to samples, as one can see 
in figures 3a and b. Both samples are almost identical 
with approximately equal mobility. Note that the Fermi 
level still lies in the bulk gap and transport is dominated 
by edge states, because we see nonlocal effects. In the 
hole side of the peak we do not find neither small neg¬ 
ative magnetoresistance near B=0, nor the oscillations 
on the side branches of the the positive magnetoresis¬ 
tance. In magnetic fields above 3 T the magnetoresis¬ 
tance falls off rapidly marking a pronounced crossover to 
the quantum Hall effect regime in accordance with pre¬ 
vious observations^^. 

The figure 4 shows the low-fled part of the relative 
magnetoresistance for the two values of the gate voltages 
and two temperatures. One can see that for this particu¬ 
lar sample the MR varies linearly with magnetic field. As 
voltage Vg increases, the B-slope of the MR decreases and 
additional oscillation emerges. The MR profile does not 
show any significant temperature dependence. Figure 4 b 
displays traces of the resistance in nonlocal configuration 
for device B {Rnl = Ri= 2 , 6 -,v= 3 , 5 )- One can see simi¬ 
lar triangular-shaped MR peak, as in the local geometry, 
though less wide, with two satellite peaks. Coexistence 
of the low-field negative MR peak in nonlocal configu¬ 
ration exclude the possibility that this effect has a bulk 
origin. It is worth noting , however, that around B=0 
the magnetoresistance is parabolic rather than linear in 
nature. 

Finally we present the results for a device C with a 
large gate. We expect that the length of the edge states 
in this structure is determined by the perimeter of the 
sample covered by the metallic gate. Figure 5 shows 
schematic view of the sample. One can see that the cur¬ 
rent flows mostly along the edges of the side branches. 
The resistance reveals saw-tooth oscillations, shown in 
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FIG. 4: (Color online) Color online) (a) The relative local 
resistance Rl{B)/Rl{0){Rl ~ Ri=i, 4 -,v= 2 , 3 ) as a function 
of the magnetic field for the two values of the gate volt¬ 
age Vg — VcNp(V)- 0.35 (black), 0.53 (blue, red), and two 
temperatures (device B). (b) The relative nonlocal resistance 
Rnl{B)/Rnl{Q) {Rnl = i?/=2,6;V=3,5) as a function of the 
magnetic held for the two values of the gate voltage Vg — VcNP 
: 0.35 (black), 0.5 (red) (device A). Dashed lines are B-linear 
approximations. 




FIG. 5: (Color online) Color online) The local resistance R for 
device C as a function of the magnetic held for the diherent 
values of the gate voltages, (a) Vg — Vcnp (V) (from top to 
bottom): 0.1, 0.3, 0.5, 0.7. 


Fig.5, when the Fermi energy passes the electronic side 
of the resistance peak away from CNF. Note, however, 
that the evolution of the oscillations with gate occurs 


rapidly and not gradually, as in devices A and B. 


III. DISCUSSION 

We now proceed to an analysis of the data described. 
We will focus on the two models that can explain the 
magnetoresistance of one dimensional edge electrons. 
Both models considered alternative paths for the edge 
states due to bulk disorder or rough edges of the realis¬ 
tic sample. Figure 6 illustrates the helical edge states in 
a disordered 2D TI in an uniform magnetic field. The 
first scenario describes disordered spinless one dimen¬ 
sional quantum wire^. For strong enough disorder elec¬ 
tron paths frequently deviate into the bulk region en¬ 
closing AB flux before returning back to the edge, la¬ 
beled in figure 6 by number 2. The conventional WAL 
approach can be used. The average with respect to the 
different size of the loops leads to the positive linear mag¬ 
netoresistance. Second scenario relies to the localization 
of the helical edge sates due to the collective action of 
the random magnetic flux through the loops, naturally 
formed by rough edges of the sample (labeled by num¬ 
ber 1 in the Figure 6 )^. In accordance with this sce¬ 
nario immediately after magnetic field is switched on, 
edge states become localized. The magnetic field pene¬ 
trates through the random helical edge states loops (Fig¬ 
ure 6 ), which acts as magnetic flux impurity and in¬ 
troduces the backscattering between edge states on each 
boundary. Similar to the first scenario the average with 
respect to the different magnetic fluxes should be per¬ 
formed. The modeli^ predicts positive magnetoresistance 
and dependence of the inverse localization length 7 
in small magnetic field, which becomes more linear with 
increasing magnetic field. Our experimentally observed 
linear positive magnetoresistance near CNF are in good 
qualitative agreement with both models. Detailed com¬ 
parison with model^^ has been performed in our previous 
publication^i. Discrimination among two scenarios re¬ 
quires further experimental work. It is worth noting that 
both models consider ballistic transport at zero magnetic 
field, while our samples demonstrate diffusive transport. 
Despite the fact that while the both models give a sat¬ 
isfactory description of the linear positive magnetoresis¬ 
tance near CNF, the explanation of negative MR and 
AB-like oscillations in Fig.3-4 away from CNF requires a 
further elaboration of the models. 

Absence of the resistance quantization in samples with 
dimensions above a few micronsiid^i^l is another un¬ 
solved problem. One of the possible explanation is the 
fluctuations of the local insulating gap width induced 
by smooth inhomogeneities, which can be represented as 
metallic puddles or dots. The well-localized metallic re¬ 
gions along the edge have been found using scanning gate 
microscopy^ and in microwaves experiments^. Accord¬ 
ing to the model 2 i charge carrier puddles coupled to a co¬ 
herent conductor results in incoherent inelastic processes 
and modifies the ballistic transport. Therefore metallic 
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FIG. 6: (Color online) Schematics of helical edge state prop¬ 
agation along a disordered edge in 2DTI in a perpendicu¬ 
lar magnetic field. Electrons tunnel in and out quantum dot 
(puddle) created by inhomogeneous charge distribution. Elec¬ 
tron helical trajectories form loops resulted from the edge 
roughness (1), near bulk disorder (2) and metallic puddles 
(3). Bulk electron trajectories form loops inside of the pud¬ 
dle. 

puddles can lead to spin dephasing, since an electron en¬ 
tering the puddle is thermalized by dissipation and later 
on fed back into the system. Therefore ballistic coherent 
transport is expected only in the region between the pud¬ 
dles, and total 4-terminal resistance exceeds the quan¬ 
tized value. Self-averaging resistance of the sample with 
edge state dominated contribution to transport is given 
by2i: 

i (1) 

where rip is the density of the puddles, A = hv/Eg ss 
18nm is the electron penetration depth into the puddles 
(v « 5.5 X lO”^ cm/s is the electron velocity, Eg ~ 2QmeV 
is the forbidden gap), g is the dimensionless conduc¬ 
tance within the dot (puddle), d is the mean level spac¬ 
ing within the dot, L is the distance between probes 
(length of the edge states). It has been found that 
the equation 1 gives a satisfactory explanation of the 
high resistance value, obtained in experiments^^. Com¬ 
bining all parameters we calculate po = R/{^L) = 

8 X 10^ {e^ jh)!cm, which is comparable with ex¬ 

perimental value po = 15 X 10^(e^//i)/cm. The density 
of the puddles increases, and the size of the puddles is 
growing, when the Fermi level lies nearer to the con¬ 
ductive band. The Figure 6 illustrates the helical edge 
states in the presence of the all types of the disorder: edge 
roughness (1), bulk disorder (2) and metallic puddles due 
to inhomogeneous charge distributions. We ague that, 


while the edge roughness as well as the bulk disorder are 
important and play a dominant role in the region near 
CNF in magnetic field, away from CNF the role of the 
metallic puddles in MR becomes more pronounced and 
essential. WL effects in ballistic cavities indeed has been 
studied both theoretically^ and experimentally^. The 
WL peak and periodic oscillations are observed which 
were attributed to the Aharonov-Bohm effect through 
a periodic orbit within the cavities. Therefore, it is 
naturally to explain the negative magnetoresistance and 
AB-like oscillations in figures 3,4 by existence of metal¬ 
lic ballistic puddles near the edge. It is worth noting 
that the Fermi energy of electrons in puddles lies in the 
parabolic part of the energy spectrum (figure I) because 
of the low density. Therefore, one would expect that 
the spin orbit coupling is weak, and magnetoresistance 
is negative in agreement with our observations. It is 
easy to estimate the characteristic sizes of these loops 
using the period of the oscillations in fig.3-4, which is 
SB = 0.3 — 0.4 T. Then the characteristic area of the 
loops is S' = di?/$o(4>o = h/e) = 10“^°cm^ and respec¬ 
tively size is about 100 nm. Indeed this value agrees with 
estimations of the puddle size and density^ and scanning 
gate microscopy^. Note that for explanation of the neg¬ 
ative sign of the MR the role of the bulk electrons in the 
puddles is emphasized, however the bulk state and edge 
state may co-exist (figure 6). The interplay between the 
topological insulators helical states and bulk electrons 
requires further theoretical study. When the Fermi level 
moves to valence band, it is expected that the puddles 
should be occupied by the holes. Note however, that both 
dephasing and spin relaxation times for holes are found 
to be much smaller than for electrons in similar condi¬ 
tions and WAL effect should be considerably smaller in 
agreement with our observations. 

When the Fermi level lies in the conductive band, and 
transport becomes dominant by massive fermions with 
Dirac spectrum (figure 1), the weak antilocalization be¬ 
haviour has been observed^i. 

In conclusion, we observed interplay between posi¬ 
tive and negative magnetoresistance, when Fermi level 
shifts with respect to the charge neutrality point, but 
still lies inside of the gap, and transport occurs via edge 
states. We consider three contributions to the magne¬ 
toresistance: edge state penetration to the bulk, edge 
state scattering by magnetic flux formed by rough edges, 
and WL of the bulk electrons in the puddles formed by 
inhomogeneous charge distributions. The negative mag¬ 
netoresistance and AB-like oscillations are attributed to 
weak localization of the bulk electrons in the metallic 
puddles formed by fluctuation of the local insulating gap. 
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